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Abstract 

We quantize the generators of the little subgroup of the asymptotic Poincare 
group of Lorentzian four-dimensional canonical quantum gravity in the continuum. 

In particular, the resulting ADM energy operator is densely defined on an ap- 
propriate Hilbert space, symmetric and essentially self-adjoint. 

Moreover, we prove a quantum analogue of the classical positivity of energy 
theorem due to Schoen and Yau. The proof uses a certain technical restriction on the 
space of states at spatial infinity which is suggested to us given the asymptotically 
flat structure available. The theorem demonstrates that several of the speculations 
regarding the stability of the theory, recently spelled out by Smolin, are false once a 
quantum version of the pre-assumptions underlying the classical positivity of energy 
theorem is imposed in the quantum theory as well. 

The quantum symmetry algebra corresponding to the generators of the little 
group faithfully represents the classical algebra. 

1 Introduction 



Following the canonical approach to a quantum theory of gravity, a (Dirac) observable is 
by definition a self-adjoint operator on the full Hilbert space (not only on the Hilbert space 
induced from the full Hilbert space by restricting to the space of solutions of the quan- 
tum constraints) which weakly commutes with the constraint operators. Equivalently, 
an observable leaves the physical Hilbert space of solutions to the quantum constraints 
invariant. 

There are no Dirac observables known in neither classical nor quantum gravity, except 
for the asymptotically flat case where it is well-known that the Poincare generators at 
spatial infinity form a closed Dirac observable algebra. 
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In the present paper we address the question of how to quantize these generators. 
We work with the real version |lj of the originally complex connection formulation of 
general relativity . The associated real connection representation has been made a solid 
foundation of the quantum theory in the series of papers ||, f|, ||, |6|, [7], ||, [| . This rigorous 
mathematical framework which is based on the earlier pioneering work on loop variables, 
for instance [pi], |TTJ for gauge theories and gravity respectively (the latter of which was 
designed for the complex variables and thus was lacking an appropriate inner product), 
equips us with the tools necessary to ask the question of whether the various operators 
constructed are densely defined, symmetric, self-adjoint, diffeomorphism invariant and so 
forth. 

But even though this kinematical framework is available, it is still quite surprising 
that something like an ADM energy operator can actually be constructed. The reason for 
this is that in the representation under consideration the ADM energy function turns out 
ot be a rather non-polynomial function of the canonical momenta and thus it is far from 
clear how to define it. It turns out that the same technique that enabled one to define the 
Wheeler-DeWitt operator for 3+1 Lorentzian gravity ]T2L O, |T3], the Wheeler-DeWitt 



operator for 2+1 Euclidean gravity length operators [[Uj, and matter Hamiltonians 



when coupled to gravity [|T7]] can be employed to define Poincare quantum operators. 

The plan of the present paper is as follows : 

In section 2 we review the necessary mathematical background from || f|, || || |7|, ||, 



In section 3 we regularize the ADM energy operator. There are at least two natural 
orderings, one in which the operator becomes densely defined on the full physical Hilbert 
space as defined in [|Ti| and one in which it is not. Nevertheless the latter operator should 
be the physically relevant one because when restricting the Hilbert space to a subspace 
which is suggested to us given the asymptotically flat structure available, then this op- 
erator turns out to be positive-semidefinite by inspection and essentially self-adjoint on 
the physical Hilbert space. This result reveals that several of the speculations spelled out 
in and which were based on several unproved assumptions were pre-mature : after 
taking the quantum dynamics of the theory and the quantum asymptotic and regularity 
conditions on the Hilbert space appropriately into account, the quantum positivity of 
energy theorem is not violated. 

It should be said from the outset, however, that the "quantum positivity of energy theo- 
rem" that we provide rests, besides on a quite particular regularization procedure of the 
ADM energy operator which exploits the fall-off behaviour of the fields at spatial infinity 
quite crucially, on one additional technical assumption (the tangle assumption, see below) 
whose physical significance is unclear. Although it can be motivated also given the struc- 
ture available at spatial infinity, it should be stressed that without this assumption the 
positivity theorem would not hold. 

In section 4 we naturally extend the unitary representation of the diffeomorphism 
group to include asymptotic translations and rotations and compute the symmetry algebra 
between time translations and the spatial Euclidean group, that is, we verify the algebra 
of the little group of the Poincare group. We find no anomaly. As is well-known, the little 
group suffices to induce the unitary irreducible representations of the Poincare group. In 
this paper we do not, however, address the more difficult problem of how to define a boost 
quantum operator. 
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2 Preliminaries 



We begin with a compact review of the relevant notions from || [L8). The interested 
reader is urged to consult these papers and references therein. 



We assume that spacetime is of the form M = R x £ where £ has an asymptotically 
flat topology, that is, there is a compact set B C £ such that £ — B is homeomorphic 
with of R 3 with a compact ball cut out. We also assume that <9£ is homeomorphic with 
the 2-sphere. The case of more than one component of <9£ (e.g. several asymptotic ends 
or horizons etc.) can be treated in a similar way. 

Denote by a, b, c, .. spatial tensor indices and by i,j,k,.. su(2) indices. The gravita- 
tional phase space is described by a canonical pair {A l a1 E°;/n) where A l a is an SU{2) 
connection on the hypersurface £, E® is an ad(S77(2)) transforming vector density and 
k is the gravitational coupling constant. This means that the symplectic structure is 
given by {A l a (x), Ej(y)} = K5(x,y)5 b 5j. The relation with the co-triad e*, the extrin- 
sic curvature K a b and the three-metric q a b = e l a e\ is Ef = \e abc 'e^kele^ and A\ = 
+ sgn(det((e J c ))) K a i,e b j where e abc is the metric independent, completely skew tensor 
density of weight one, T l a is the spin-connection of e l a and e" is the inverse of the matrix 

By 7 we will denote in the sequel a closed, piecewise analytic graph embedded into 
a d-dimensional smooth manifold £ (the case of interest in general relativity is d = 3). 
The set of its edges will be denoted E{p/) and the set of its vertices V(j). By suitably 
subdividing edges into two halves we can assume that all of them are outgoing from a 
vertex (the remaining endpoint of the so divided edges is not a vertex of the graph because 
it is a point of analyticity). Let A be a G-connection for a compact gauge group G (the 
case of interest in general relativity is G = SU(2)). We will denote by h e (A) the holonomy 
of A along the edge e. Let 7Tj be the (once and for all fixed representant of the equivalence 
class of the set of) j-th irreducible representations of G (in general relativity j is just a 
spin quantum number) and label each edge e of 7 with a label j e . Let v be an n-valent 
vertex of 7 and let ei, .., e n be the edges incident at v. Consider the decomposition of the 
tensor product ®fc =1 7Tj e into irreducibles and denote by 7r Ctl (l) the linearly independent 
projectors onto the irreducible representations c v that appear. 

Definition 2.1 An extended spin-network state is defined by 

T ^M) '■= *K®t*=V( 7 )kc„(l) • ^eeE^vee^MA))]) (2.1) 

where j = {j e }e€E(j), c = {cv} v eVM- I n what follows we will use a compound label 
I = (j(I),j(I),c(I)). An ordinary spin-network state is an extended one with all vertex 
projectors corresponding to singlets. 

Thus, a spin-network state is a particular function of smooth connections restricted to a 
graph with definite transformation properties under gauge transformations at the vertices. 
Their importance is that they form an orthonormal basis for a Hilbert space 7i = 7i aux , 
called the auxiliary Hilbert space. Orthonormality means that 

< -^yj,c> -^y'J',? > = < ^yj,?' -^y'J',? > aux= Sjry'SjjfSg^ . (2.2) 
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Another way to describe Tt is by displaying it as a space of square integrable functions 
L2(A/G,dp,o). Here A/Q is a space of distributional connections modulo gauge transfor- 
mations, typically non-smooth and fio is a rigorously defined, a-additive, diffeomorphism 
invariant probability measure on A/Q. The space A/Q is the maximal extension of the 
space A/Q of smooth connections such that (the Gel'fand transform of) spin-network 
functions are still continuous. The inner product can be extended, with the same or- 
thonormality relations, to any smooth (rather than analytic) graph with a finite number 
of edges and to non-gauge-invariant functions. It is only the latter description of TC which 
enables one to verify that the inner product < ., . > is the unique one that incorporates 



the correct reality conditions that A, E are in fact real valued. The inner product (2.2) 
was postulated earlier (see remarks in [p9fl ) for the complex connection formulation. But 
it was not until the construction of the Ashtekar-Lewandowski measure /io that one could 
show that this inner product is actually the correct one for the real connection formulation 
only. 

We will denote by $ the finite linear combinations of spin-network functions and call 
it the space of cylindrical functions. A function / 7 is said to be cylindrical with respect 
to a graph 7 whenever it is a linear combination of spin-network functions on that graph 
such that TTj e is not the trivial representation for no e G Ei^y). The space $ can be 
equipped with one of the standard nuclear topologies induced from G n because on each 
graph 7 every cylindrical function / 7 becomes a function /„ on G n where n is the number 
of edges e of 7 through the simple relation f^(A) = f n (h ei (A), .., h en (A)). This turns 
it into a topological vector space. By $' we mean the topological dual of $, that is, 
the bounded linear functionals on $. General theorems on nuclear spaces show that the 
inclusion $ C H C $' (Gel'fand triple) holds. 

So far we have dealt with solutions to the Gauss constraint only, that is, we have 
explicitly solved it by dealing with gauge invariant functions only. We now turn to the 
solutions to the diffeomorphism constraint (we follow [|H|). 

Roughly speaking one constructs a certain subspace $d«// of $' by "averaging spin- 
network states over the diffeomorphism group" by following the subsequent recipe : 
Take a spin- network state Tj and consider its orbit {T/} under the diffeomorphism group. 
Here we mean orbit under asymptotically identity diffeomorphisms only ! Then construct 
the distribution 

[7}]:= £ T (2.3) 
re{Tj} 

which can be explicitly shown to be an element of Any other vector is averaged by 
first decomposing it into spin-network states and then averaging those spin-network states 
separately. Certain technical difficulties having to do with superselection rules and graph 
symmetries || were removed in 



An inner product on the space of the so constructed states is given by 

<[fU9]>Dif f :=[f](9) (2.4) 

where the brackets stand for the averaging process and the right hand side means eval- 
uation of a distribution on a test function. The completion of $Di// with respect to 
< ., . >Diff is denoted H D iff- 



Finally, the Hamiltonian constraint is solved as follows [14| : One can explicitly write 



down an algorithm of how to construct the most general solution. It turns out that one 
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can construct "basic" solutions G $' which are mutually orthonormal with respect to 
< .,. >Diff (in a generalized sense) and diffeomorphism invariant. The span of these 
solutions is equipped with the natural orthonormal basis s M (in the generalized sense). 
One now defines a "projector" 

^/:=[[/]]:=E a M<*M. [/]>»// ( 2 - 5 ) 

for each / e <£> and so obtains a subspace <&Ham C The physical inner product is 
defined by 

<[[/]UM]>^:=[[/]]([9]). (2-6) 

Finally, the physical Hilbert space is just the completion of $# am with respect to < 

• j • ^ Ham- 



3 Regular izat ion of the ADM Hamiltonian 

There are many ways to write the ADM-Hamiltonian which are all classically weakly 
identical. We are going to choose a form which is a pure surface integral and which 
depends only on Ef because in this case the associated operator will be almost diagonal 
in a spin- network basis so that we can claim that spin-network states really do provide a 
non-linear Fock representation for quantum general relativity as announced in O, jl3fl. 



The appropriate form of the classical symmetry generators was derived in |2(| |21 
Although that paper was written for the complex Ashtekar variables, all results can be 
taken over by carefully removing factors of % at various places. We find for the surface 



part of the Hamiltonian (expression (4.31) in |2(J, we use that N = iV/ydet(g), D a N = 
(D a N) / det(q) where N is the scalar lapse function) 

Or AT 

E(N) = — / dS a . Efd b E b . (3.1) 
kJ 9x v /det(g) 

It is easy, instructive and for the sign of the ADM energy crucial to see that ( |3.l| ) really 
equals the classical expression + - Jg-^dS a (q abjb — Qbb,a) due to ADM : Using that Ef = 
^e abc 6ijke J a e% we have the chain of identities 

Efd b E b = -sgn(det(e))e^% fc [e^], 6 



det(g) 



-2sgn(det(e)Ke 6c %^< 6 

-2sgn(det(e))g a ^% fc e}e^ = -2q a f e bcd ^^)e fce ele% 
-Aq a %5 d e] ^t{q)etel b 

A^^)q ac q e %el e] = 2^t{qjq ac q b %[el b - eQ 
^^)q ac q bd [2e\ d e^ b + 2e]/^ - 2e\ d e\ U ~ HAA 
det(q)q ac q bd [(q cd , b - q bd , c ) + 2e\ d e\ b ] (3.2) 
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Now we expand e l a {x) = S l a + ^^J^ + o(l/r 2 ) where r 2 = 5 ab x a x h defines the asymptotic 
Cartesian frame. The function f l a {x/r) only depends on the angular coordinates of the 
asymptotic sphere and is related to the analogous expansion q ab (x) = 5 ab + ^ ab ^ r ^ +o(l/r 2 ) 
by f ab 5i = /*. Consider now the remainder in the last line of (|3.2|) . Its o(l/r 2 ) part 
vanishes because /[ a &] = and this concludes the proof. 

In the sequel we focus on the energy functional Eadm = E(N = 1). We will quantize 
it in two different ways correspoinding to two quite different factor orderings. Each of the 
orderings has certain advantages and certain disadvantages which we will point out. 

3.1 Ordering I : No state space restriction 

In this subsection we will derive a form of the operator which is densely defined on the 
whole Hilbert space TC (and extends to the spaces Hniff, 'Hphys defined above) without 
imposing any further restriction that corresponds to asymptotic flatness. 
Using again that E^ = \eijk^ abc e b e k we can write it as 

Eadm = lim E ADM (S) where E ADM (S) = -- [ 1 e ij V A e k d h E\ (3.3) 
s-es K Js ^det(g) 

and S is a closed 2-surface which is topologically a sphere. The idea is to point-split 
expression (|3.3| ) and to use that 

[sgn(det(e))4](x) = ±{j£(x), V(x, e)} 



where V(x,e) = J E d 3 yXe(x,y)ydet(q)(y) and Xe is t ne (smoothed out) characteristic 
function of a box of coordinate volume e 3 . Since 

lim ^ ' ^ = 5(x, y) so that lim — - J - = \/det(q)(x) 



we have 



Eadm(S) 

lim-- / . 1 ^V(x)Ae fe (x) / d 3 y X e(x, y){d h E*){y) 



^ o k Js V [x, e) js 

lim n -^~3 I 7^r—AA 3 (x),V(x,e)} A{A k (x),V(x,e)} [ d 3 y X e(x,y)(d b E b )(y) 
^ o 2k 6 Js V{x, e) js 

lim — / <# k {Ai(x), y/V(x, e)} A {A k (x), y/v&, e)} / dVfcfo y)(ftJS?)(y) 

e— J5 v v JS 

liml | tr({A(x), ^(x, e)} A {A(x), ^V(x, e)} J^d 3 yXe(x,y)(d b E b )(y)) 
-lim 1 / tr({A(x), ^(s, e)} A {A(a;), ^(a, e)} / d 3 y[d ybXe (x, y)]E b (y)) 

e->0 K Js v v JS 

limi% DM (£) (3.4) 
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where in the second before the last step we have taken a trace with respect to generators Tj 
of su{2) obeying [n, Tj] = e^Tfe and in the last step we have performed an integration by 
parts (the boundary term at <9E does not contribute for finite 5* and e sufficiently small). 



Thus, we absorbed the l/ydet(g) into a square- root within a Poisson-bracket and simul- 
tanously the singular 1/e 3 into a volume functional. Classically we could have dropped 



the l/ydet(g) (although the integrand would then no longer be a density of weight one 
and is strictly speaking not the boundary integral of a variation of the Hamiltonian con- 
straint) due to the classical boundary conditions which tell us that det(g) tends to 1. 
We now quantize E ADM (S). This consists of two parts : In the first we focus on the 



volume integral in (3^) and replace Ef by Ef = —iknS/SA^. In the second step we trian- 
gulate S exactly as the hypersurface of 2+1 gravity in [|T5[], replace the volume functional 
by the volume operator and Poisson brackets by commutators times l/(ih). 

So let / 7 be a function cylindrical with respect to a graph 7. Since we are only 
interested in the limit S — > <9£ we may assume that 

1) 7 lies entirely within the closed ball whose boundary is S but 

2) 7 may intersect S at an endpoint of one of its edges and may even have edges that lie 
entirely inside S. 

Furthermore we can label the edges of 7 in such a way that an edge either intersects S 
transversally (with an orientation outgoing from the intersection point with S) or lies 
entirely within S. 

Coming to the first step we have for the y integral involved in E ADM (S) : 
d z y[d a x,{x,y)}E a ( y) i f 1 
= -ihn I d 3 y[d aX e(x,y)} [' dte a (t)S(y,e(t))Xi(t)f 7 

= -ihK / dte a (t)[d y aXe(x,y) y=<t )]Xl( t )fj 

ee£( 7 ) ° 

= - ihK £ f l dt[^-Xe(x,e(t))}Xi(t)f 7 
eeS ( 7 ) J o at 

n 

= -ihK ti^Z!ta( x > e (*fc)) ~ Xe(z> e(*fc_i))]X*(*fc_i)/ 7 (3.5) 

e£E(-y) ^ k=l 

where E(j) is the set of edges of 7, X l e {t) := [7i e (0, t)rih e (t, l)]ABd/d[h e (0, 1)]ab and 
= to < ti < •• < tn — 1 is an arbitrary partition of the interval [0, 1]. 
It is important for what follows that for each t,e,i X*(t)/ 7 is still a function cylindrical 
with respect to 7. 

We now come to the second step. This involves, first of all, a triangulation of the 
two-dimensional surface S in adaption to the graph 7. Besides the prescription explained 
in detail in |15| which deals with the triangulation of S in the neighbourhood of a vertex 



formed by edges of 7 that lie entirely within S we just need to deal with the case that 
a vertex v of 7 has also edges incident at it which lie entirely inside the open ball with 
boundary S except for the one point v. In case that there are at least two edges ei,e 2 



of 7 incident at v such that ex, t<i C S we can still take over the triangulation from |L5 . 
However, if there is only one or no such edge (indeed, since we do not allow for gauge 
transformations at spatial infinity we can allow for open edges that lie entirely within or 
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end at S without ruining gauge invariance) we need an additional prescription : In case 
there is only one edge e C S incident at v, choose an arbitrary edge e! not intersecting 7 
except at v such that the tangents of the edges e, e', e" are positively oriented at v where 
e" is any of the edges of 7 incident at v but transversal to S. 

In case there is no edge e C S incident at v, choose two arbitrary edges e, e' not intersecting 
7 except at v such that the tangents of the edges s, e', e" are positively oriented at v where 
e" is any of the edges of 7 incident at v but transversal to S. These arbitrary edges will 
disappear from the stage again at the end of the calculation. 

Given this set-up, at each vertex v of 7 that lies inside S we have now at least to edges 
ei,ea incident at it that are inside S and we can define the triangles A associated with 
pairs of edges incident at v and inside S exactly as in As in |15|] we then have 



two segments Si(A), 52(A) for each triangle A which are actually segments of edges of 7 
incident at v and that lie inside S. Now observe that 



^h Si (A){h s \ A) , V(v, e)}h s . {A) {h Sj \ A) , V(v, e)} 



= 5V^0)^(0){,4>), y/V(v,e)}{A(v), y/V(v,e)} + o(6 2 ) 

= 2vo\(A)e ab {A a (v), ^Jv(v, e)}{A b (v), ^Jv(v, e)} (3.6) 

where 5 is a small parameter corresponding to the parameter length of the Sj(A) and 
e ab the metric independent totally skew tensor density of weight one on S. Altogether 
we therefore conclude that the surface integral J s {A(x), \jv(x, e)} A {A(x), ^Jv(x, e)} 
involved in Eadm(S) can be quantized by 

E £73 E ^h aiW [h~\ A) , y/v(v,e)]h 8j(A) [h-\ A) , y/v(v,*)]fy ■ (3-7) 

«gV(7) v ' v(A)=v 



We do not wish to give a full derivation of ( |3.7| ) for which the reader should consult fl5 | 
however, a few remarks are in order which intuitively explain 



1) First of all as already mentioned, we do not work with a fixed triangulation but 
with a whole finite family of triangulations that depend on the graph 7 of the state 
that we act on. More precisely, a particular member of the family of triangulations 
is defined such that for each vertex v of 7 which lies in S we a) pick one pair ei, 
of edges of 7 incident at it which are such that ti is next to t\ to its right while the 
tangents of ei,e2 enclose an angle less than or equal to n (with respect to 5 a b), b) 
take proper subsegments Sj(A) C incident at v, c) construct two more segments 
si(t) = 2v — Si(t), d) construct four obvious triangles from Sj, which saturate v and 
e) choose a triangulation which embeds those four triangles for each v as basic ones 
and is otherwise arbitrary only subject to the restriction that none of the remaining 
A has its basepoint on 7. 

2) Next, each of the four triangles Ai 2 , A 2 i, Ajg, Agi contributes classically the same 
to the surface integral f s — a I a as ^ ne triangle Ai 2 so that we can classically 
replace these four terms by 4 J Al2 - This explains the factor of 4 in ( |3.7| ). 

3) Now, since we do not want to distinguish one particular pair of edges as compared 
to any other, we average over the choice of pairs which explains the factor of 1/E(v) 
where E(v) = n(v) or E{v) = n(v) — 1 is the possible number of pairs depending 
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on whether the angles less than or equal to 7r add up to 2% (with respect to 8 a b) or 
not and n(v) is the valence of v. 

4) Finally, the fact that we sum over vertices of 7 only comes from the presence of 
the volume operator which has non-trivial action at vertices only. Therefore the 
contribution of all other triangles which define the triangulation drop out (all tri- 
angles such that there is no edge of 7 transversal to 5* at its basepoint since the 
three-dimensional volume operator annihilates co-planar vertices). 

It should be mentioned that any of the so-defined "averaged family of triangulations of S 
adapted to 7" has as classical continuum limit the original integral J E ! 
Putting Q3.5p and ( |3.7| ) together we obtain as the final result 

v£V{y) v ' v(/\)=v 

x e ij tv(h Si(A) [h~ l (A) , \Jv(v,e)]h SjiA) [h~\ A) , \Jv(v,e)\ x 

x E i2\Xe(v,e(t k ))- Xe (v,e(t k ^))}Xl(t k ^)f,. (3.8) 

e£E(-y) k=l 

Now we perform the limit n — > 00 and e — > in reversed orderf] : Keeping n fixed, for small 
enough e only the term with k — 1 in the sum survives provided that e(0) = v. Therefore 
also the sum over k and with it the n dependence drops out. Also the operator V(v, e) 
actually has a limit as e — > which we call V v and which is defined by lim^o V(v, e)/ 7 =: 
(14)7/7 f° r au 7- That the limit exists relies on the fact that either 7 has v as a vertex 
or it does not. In the latter case the limit just vanishes, in the former case for sufficiently 
small e the e box around v does not include any other vertex of 7 other than v and so 
there is only one contribution (Vu) 7 / 7 which is constant as e —>■ 0. The family of operators 
{y v )~f is consistently defined because V is. For a more explicit formula in terms of analytic 
germs of edges see [[TJJ. 

Now, the limit n — > 00 is trivial and the resulting operator derived for arbitrary but finite 
S can be extended to <9£. The result is 

EADMfn, = jr^ E 7^ £ x 



x e ij tr(/i Si(A) [/i s .; A) , >/^e)m(A)[/i,/ ( A). V^>^)] x 

E (3-9) 

eeB(7),e(0)=« 

where X l e (0) = X l e = X % {h e )) and X i (g) is the right invariant vector field at g G SU(2). 
The virtue of ( p.9|) is that it displays the ADM energy operator as a densely defined 



1 One can also make e y-dependent in ( |3.4| ) in a state-dependent way which then leads to a dependence 
e = o{l/n). If one then takes n — > 00 one gets classically back to the original expression for the ADM 
energy. If one takes n large but finite then one arrives at the same result as below on the quantum level 
without that an interchange of limits is necessary. Therefore the calculations that follow are justified. The 
state-dependence of the regularization drops out in the final expression because before actually taking 
the limit n->oo the operator will be n-independent. The family of operators so obtained for each state 
are consistently defined as we will see. See @ for a more detailed explanation. 



9 



operator on all of the Hilbert space. Also, the dependence on the "arbitrarily short edges 
Si (A) drops out at the end of the calculation because of gauge invariance as explained in 
fl6|j . The disadvantage is that the operator ( |3.9| ) is not a manifestly positive semi-definite 
operator. This is, however, not surprising because even the classical ADM energy is not 
a positive semi-definite functional on the full phase space of general relativity. It is only 
when evaluating it on a) asymptotically flat b) solutions of the Einstein equations which 
c) satisfy an energy condition for allowed matter and d) allow for a regular initial data 
set, that the positive energy theorem has been proved |22| , |23|| and in fact one can easily 
produce negative ADM energy when one of these conditions is violated. As we did not 
impose any (quantum analogue of) such restrictions we cannot expect to find a manifestly 
non-negative operator. 

In the next subsection we will derive another quantization of the ADM energy which is 
only densely defined on a subspace of the Hilbert space, however, the definition of that 
subspace a quite natural quantum translation of the classical condition that there be an 
asymptotically flat regular initial data set. The virtue will be that the ADM operator 
acquires non-negative discrete spectrum on that subspace of the Hilbert space thus proving 
a "Quantum Positivity of Energy Theorem" . 



3.2 Ordering II : Restrictions on the State Space 

In order to make sense of the operator to be defined in this section we need to give 
some definition of "asymptotically flat state". The following definition is a first attempt 
towards a precise notion of "quantum asymptotic flatness" to be considerably refined in 
future publications. 

Definition 3.1 An asymptotically flat state \I/ on the Hilbert space 7i is a distribution in 
$' satisfying the following conditions : 

i) ^ is a normalized solution fj^ to all quantum constraints U(<p) — 1 = 0, H(N) = (the 
diffeomorphism and Hamiltonian constraint) of general relativity where if, N are arbitrary 
three- diffeomorphisms and lapse functions subject to the condition that they be pure gauge, 
that is, <p(x) — > id and N(x) — > as x — > <9£ and where fj is the operator that maps 
elements of $ to solutions to all constraints. 

ii) ^ — fj^ is asymptotically flat. That is, consider any compact region R, surface S 
and loop a ( which is embedded in the graph underlying the definition of of £ — <9S 
and which are large as compared to Planck volume, area and length respectively as mea- 
sured by the Euclidean metric := 5 a b- Then, as R,S,a tend to <9£ ; the quanti- 
ties V([V(R) - V Q (R)]Z),V([A(S) - A (S)}0,y([tr(h a ) - 2]f) respectively are of order 
£p, £p, tp/ L(c) respectively where Vo(R) , A (S) , L Q (a) are the volume, area and length of 

R, S, a as measured by q^ . 

Hi) \l> transfroms according to an unitary, irreducible representation of the Poincare group 
at spatial infinity. 

iv) ^ satisfies the dominant energy condition : Let H' matter (N) be the dual of the matter 
part of the Hamiltonian constraint (not Hamiltonian !) and V^^^N) be the dual of the 
matter part of the Diffeomorphism constraint. 

The four-vector = (N, N) is said to be a future directed timelike vector in a state ^> 
if a) N > and b) there exists e > such that -t 2 N 2 (x) + ty[L(c{N,x,t)) 2 £] < for 



each < t < e where L(c) is the length operator [TqJ and c(N,x,t) is the segment of the 
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integral curve of N beginning in x and ending after a parameter distance t. 
\I/ is said to satisfy the dominant energy condition provided that for every four-vector 
which is future directed and timelike for \l/ then there is an e > such that the four vector 
P^(x,t) := (H' matter (N x j) + V^aiier(-^,t)) 0) ^ s either zero or future directed and timelike 
for every x£E and < t < e in the state \l/ where N x j(y) = Xt{ x i y)N{y) and likewise for 
N x ,t{y) (Xt{x, y) is the characteristic function of a box of coordinate volume t 3 and centre 
x). In other words, ^([Hmatter{N X;t ) + V maUer (N X)t ))C) > for each x, < t < e. Here we 
have adopted the convention that the signature of the Lorentz metric be —,+,+,+. 
The subspace of £ 's inTt satisfying these conditions will be called 7i a f where u af" stands 
for asymptotically flat. 

Condition i) makes sure that \I/ is a solution of the "Quantum Einstein Equations" . Con- 
dition ii) is a possible way of defining asymptotic flatness q a b — > Sab, Kab — ► (although 
not very carefully, no fall-off and parity conditions were imposed |20], [21]] and certainly this 



condition needs to be refined in future publications. For instance, one might imagine that 
the error of ^ (V (R)£) /Vq(R) — 1 is even smaller than £p/Vo(R) in the sense that it could 
depend on some negative power of the value of the radius r (with respect to 5 a b) at the 
center of R). Condition iii) makes contact with physics and allows us to identify certain 
states with elementary particles. In particular, in the present context of pure gravity we 
should be able to isolate the graviton (spin 2, massless) states. Notice that in order to 
allow for non-trivial representations of the little subgroup of the Lorentz group we must 
specify the appropriate Diffeomorphism group in the group averaging process in order 
to arrive at the diffeomorphism invariant states |J which means, roughly speaking, that 
we include only those diffeomorphisms in the averaging process that approach identity at 
<9£. Finally, Condition iv) is the most imprecise one and attemts at defining a possible 
quantum analogue for the dominant energy condition : recall that the classical energy 
momentum tensor is said to satisfy the dominant energy condition if for every future 
directed timelike vector field i> M , the vector field T^v v is zero or future directed and time- 
like. Now if V 1 = Nn 11 + iV M is the future directed timelike foliation vector field underlying 
the split M = RxE and n M the normal vector field of £ then we may pick a frame such 
that n tJ- T t j -u t u = NH maU + N a {y ma tt)a is the only non- vanishing component of the vector 
T^ v t v and so we need to ask that it be non-negative, in particular, for N = 0, we ask 
that the matter Hamitonian densities be non-negative. This condition is, of course, in- 
complete since it is frame dependent and needs to be improved in the future. Notice that 
ifg^tn u )(x) = -t 2 N 2 (x)+t 2 (q ab N a N b )(x) w -N 2 (x)t 2 + L(c(N, x, t)) 2 which motivates 
our definition of future directedness of t M . We conclude with the remark that Condition 
iv) is certainly satisfied in the vaccum case that we are interested in in the present paper. 



One might think that a state satisfying all those conditions is rather hard to construct. 
Let us pause for a moment to argue that it is rather simple : 

Consider for simplicity a S with topology of R 3 and distribute a countable number of 
vertices v n randomly into £ with an average next neighbour distance of a Planck length as 
measured by q a ° b . Make v n the only vertex of a graph 7„ which is four-valent and non-co- 
planar (for instance, 7„ = a n fl (3 n where a n , f3 n are two kinks with vertex v n ). The graphs 
7„ are supposed to be contained in a box B n of volume ki 3 for some positive number 
k and the B n are mutually non-inersecting. Consider normalized vectors /„ which are 
finite linear combination of spin-network states defined on 7„ and which are eigenstates of 
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the volume operator V(R) for any region R, all with the same eigenvalue X n £^ = X£p > 
if v n G R. Thus we have V(B m )f n = X£p8 min f n . Consider the infinite product state 
£ := IT^Li fn which is a regular (non-cylindrical) spin- network state on the infinite graph 
7 = U n 7„ and which is in fact normalized, 1 1£| | = 1 thanks to the disjointness of the graphs 
7 n because of which ||£|| = fin ||/n|| due to the properties of the Ashtekar-Lewandowski 
measure. We now choose k := A and find that for any macroscopic R, that is, any R 
that contains many of the boxes B n it holds that V(R)£ = Vo(i?)[l + o(£p/Vo(R)]£. Now, 
since no state which is cylindrical with respect to any of the graphs 7„ can be in the 



image of the Hamiltonian constraint [12, 13, fj it follows from its definition [18 that the 
rj operator reduces to group averaging with respect to the diffeomorphism group because 
of which the group averaged diffeomorphism invariant state \I/ = [£] = r/£ is normalized 



as well with respect to the physical inner product |L8j 11^11?^ = ^(0 = ||£| 



Thus indeed V([V(R) - V (R)]£) = o(£%/V (R)) is satisfied. It is clear that the con- 
struction can be repeated for the surface operator as well because most of the inter- 
sections of the macroscopic surface S with the 7„ will not be in vertices of the 7„ so 
that V(R), A(S) can be simultanously diagonalized up to errors of order of £p/A (S). 
Thus, almost every /„ can be chosen as a simultanous eigenvector of V(R),V(S). Fi- 
nally, any macroscopic, for simplicity non-self-intersecting (any loop is product of these), 
loop a on our particular 7 is of the product form a = ° n a^ n , a n C j n , k n G {0,1} 
where k n = except for finitely many. The SU(2) Mandelstam algebra is too com- 
plicated as to exhibit an explicit solution for SU(2) so let us argue with an U(l) sub- 
stitute that the condition stated in definition (|3.1| ) is reasonable. For U(\) we have 
h a = Ilfc n =i h an - Now, if we choose for simplicity a n = j n then f n = J2k=-N a kha n where 
Xk(g) = 9 k is the character of the irreducible representation of U(l) with weight k. Since 
T = XkXi — Xk+i the condition stated in the definition amounts to asking that (for U(l)) 
1 = llfc n =i Y,k \ a k\ 2 = rifc^i Y.k=-N+i ^fc^fc+i U P to some corrections. Indeed, if we could 
choose all au to be equal (= l/\ / 2N + 1) then the error would be 1 — rifc n =i[l — 1/ (2iV+l)] 
which is small provided that Y,k n =i 1 = (.L{a)/l p ) « N. For instance we may choose 
iV = [L/£ p ] 2 where L is the bound on the length of a macrosccopic loop that we wish 
to consider. This will suffice to motivate definition ( |3.1| ). Obviously one has to refine it 
but this seems impossible without the notion of coherent states and will be left for future 
publications 



We will now show that on an asymptotically flat state the ADM energy operator as 
defined below is non-negative. We will not show that vanishing energy corresponds to 
Minkowski space. Our definition of asymptotically flat states as of yet is not restrictive 
or precise enough for that. As we will see, in order for the ADM operator to be densely 
defined we need that the following stronger condition : 

ii) 7 A state ^ = fj^ is said to be asymptotically flat provided ii) of Definition \3. i| holds 
and in addition : 

Let 7 be the (infinite) graph on which £ depends, p a point in <9£ fl ^(7) and B t , t G [0, 1] 
any homotopy of regions in S such that p G B t for each t and B = {p} . Then we require 
that for each such p there exists e > such that £ is, for each < t < e, a finite linear 
combination of eigenstates with non-vanishing and t— independent eigenvalues of the vol- 
ume operator V (Bt) . 
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It is not clear that condition ii) implies ii)' and if that should not be the case then 
we must add the requirement ii)' stated as an additional restriction on H a f ! Notice, 
however, that ii)' is not unreasonable in the asymptotically flat context. 
It will turn out in the course of the derivation that the positivity of energy theorem then 
holds if we impose one additional condition on the so already restricted space of states. 

We write expression ( |3.4| ) this time in the form (setting S = <9£ right from the begin- 
ning) 

E ADM = lim--/ dS a Ef(x) = - / d 3 y Xe (x,y)(d b E*)(y) 
= lim-- / dS a E?{x)-^— I d 3 y X e(x,y)(d b E^(y) 

t^O K JOT, V [x, e) JT, 

= lim-- f dS a (x)Et(x)—^— f d s y X e(x,y)[Gi(y)-e ijk Ai(y)E b k (y)} 

t^O k JdT, V [x, e) JT, 

=: limEW (3.10) 

e— hj 

where Gi = d a E£+ [A a , E a ]i is the Gauss law constraint. Recall that G{ = only needs to 
hold in the interior of £ because the Lagrange multiplier A* of the Gauss constraint falls 
off like 1/r 2 so that at <9E every function of Ef, A l a is gauge invariant. More precisely we 
have the following : It is of interest by itself to derive the quantum Gauss law operator 
on a function of smooth connections cylindrical with respect to a graph 7 : 

d 3 xA i (x)G i (x)f 7 



%k I ' d z xA\x) f ^^e a (t)([^^ + e^ a (x)]5(x,e(t)))(X e fe (t)/ 7 ) 

ee£(7) ° 

%k £ f d 3 xA\x) f 1 dt([-^6 ik + e tjk e a (t)Ai(x)]6(x,e(t)))(^(t)f,) 



+ f dtS{x,e{t)){Xl{t) + e ljk e a (t)Ai(x)X k e (t)}f J } 

J 

= -i%K{ £ [-A*(e(l))X:(l)+A*(e(0))X:(0) 

eS-B(7) 

+ f dtA\e(t))(Xl(t) + 6 ufc e a (t)^ a (e(t))X e fc (t)]/ 7 } . (3.11) 

<J 

Here we have made use of d x d~(x,y) = —d y S(x,y) in the third step which holds on spaces 
of test functions of rapid decrease. But since for smooth connections h e (t,t + St)^ 1 = 
1 ± e a {t)StA a (e(t)) + o(5t 2 ) we have 

X l e (t + St) - Kit) = tr(/i e (0, t + St)nh e (t + St, l)d hem ) - X\t) 
= tr(/i e (0, t)h e (t, t + St)nK% t + 5t)h e (t, l)4 e( o,i)) - x\t) 
= ^e a (t)tr(/i e (0,t)[A a (e(0),Ti]/i e (t,l)4 e(0ll )) + o(^ 2 ) 
= -e i3k 5te a (t)Ai(e(t))tr(h e (0, t)r k h e (t, l)d hei0 ,i)) + o(5t 2 ) 
= -e l3k 5te a {t)A{{e{t))X k e {t) + o(5t 2 ) . (3.12) 
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This shows that the t— integral in ( |3.11| ) vanishes identically for smooth connections. 
Now Xl(0) = X l {h e ) =: X\ where X i (g) is the right invariant vector field at g G SU(2) 
and X l e (l) = — X % (h~ l ). Thus, when splitting each edge into two halves e = e\ oe^ 1 where 
both ei,e2 are outgoing at the vertex different from their intersection point then e(l) = 
e2(0) and the right invariance of X now implies X(h e ) = X(h ei ), —X(h~) = X(h e2 ). 
Summarizing we find that the Quantum Gauss Constraint is given by 

6(A)/ 7 = -z/we £ A*(e(0))Xg/ 7 = —ifiK £ A i (v)X i v f 1 (3.13) 

where — iX v = —i J2 e (o)=v X e is the total "internal" angular momentum operator. Notice 
that ( |3.13j ) can be extended from smooth to distributional connections and that no as- 
sumption on asymptotic behaviour or smoothness of A had to be made. The Quantum 
Gauss Constraint is obviously a self-adjoint operator on 7i and anomaly free : it is trivial 
to check that [G(A),G(A')] = G([A,A']) precisely mirroring the classical constraint alge- 
bra. Moreover, the Quantum Gauss law constraint is identically satisfied as 7 tends to 
<9E because A^s = 0, that is, there are no internal charges in general relativity 0. This 
allows quantum states of distributional connections to be non-gauge invariant at spatial 
infinity, a fact that we are going to exploit in the sequel. 

Let now the state £ G Ti a f be considered as a function / 7 cylindrical with respect to a 
graph 7 which is a finite subgraph of the graph on which £ depends and which inersects 
<9£. Because £ G H a f we know that / 7 is a finite linear combination of eigenstates of the 
volume operator V(R) with non-zero eigenvalue for sufficiently small regions R and such 
that i?nl/(7)|as 7^ 0. Consider first the volume integral in fl3.10|) . Setting A(y) = Xe( x , y) 



there is an obvious quantization for the term proportional to Gi in view of ( |3.13|) . However, 
for the remainder we have, again on functions of smooth connections only to begin with 



£ijk y s d s y Xe (x, y)Al(y)E%(y)f 1 

= -ihKe ijk £ [ d s y X e(x,y)Ai(y) [* dte a (t)5(y,e(t))X*(t)f 7 
ee£( 7 ) 

= -ihKe ijk £ f 1 dte a (t)Ai(e(t)) X e(x,e(t))X k e (t)f, 
ees( 7 ) Jo 

n 

= 2ihn J2 J™ E Xe^, e^-^tr^f/ie^-i, it) - 1]) jq(t*-i)/ 7 . (3.14) 

ee£(7) k=l 

Now recall that x G <9£ and that 7 C B(dT>). Pick a particular edge e in fl3.14j ). Then 
for sufficiently small e the corresponding term either vanishes or e(0) = x and the term 
becomes tx(ri[h e (0, t(e)) — l])X*/ 7 + o(e 2 ) where t(e) is the largest value of t such that 
Xe{x,e(t)) = 1. Now for a classical, smooth connection which approaches infinity as 1/r 2 
the term h e (0,t(e)) — 1 is at most of order e/r (the "length" of e(e) is at most of order 
r) and so vanishes even at finite e because x — > <9E. Therefore expression ( |3.14j ) vanishes 
and the volume integral contribution of Eadm becomes 

/ d 3 y X e(x,y)d a E?(y)f, = -thK £ X e(x,v)X^ (3.15) 
which one can extend to non-smooth connections. 
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We turn to the surface integral of ( |3.10| ) and write f 7)V = X*/ 7 We have, ordering the 
1/V(x,e) to the left 



/ Xe(x,v)dS a y 1 - Eff J>v 
J9T, v(x,e) 



ihK [ L dte a (t)[ dS a (x)S(x,e(t))^^-Xi{t)f J>v 

*r; x Jo JdT, V x. f) 



ee£(7) 



as ' ' V(x,e) 



= —inn sgn(aL,e) — X e j 1>v 

V(e(0,e " 

where sgn(S', e) is the sign of the intersection of e with the surface S (which is outward 
oriented) at e(0) which is thus —1 because all edges e are outgoing from a vertex e(0) 
and, because of 7 C i?(<9£), they are thus running away from <9E (there is no contribution 
from edges that run inside <9E because sgn(9E, e)\ e (t) — for all t as was shown in [[25] ). 
Thus, only edges which run transversally into contribute to the sum in ( |3.16| ). 
We can take now the limit e — > 0. Notice that for small enough e we have Xe( e (0),f) = 1 
as e — > provided that e(0) = v is a vertex of 7. Moreover, for small enough e, t> is the 
only vertex of 7 in the e— box around v. Thus, we may replace V(v, e) by V v interpreting 
the operator 1/V V by its spectral resolution. Now the only critical point is whether the 
two operators X l e1 X l v that are to the right of l/V v will leave the crucial property of / 7 
intact, namely, that / 7 = fi with V v fi = Xifi, Aj 7^ 0. But this is easily seen to be 
the case provided we impose the following additional restriction on TC a f : 

Tangle Property 

An asymptotically flat state is a linear combination of diffeomorphism group averaged el- 
ements of the kinematical Hilbert space Ti a f each of which depends on an (infinite) graph 
7 all of whose intersections with <9£ are transversal, that is, there are no edges 0/7 which 
lie inside <9£. Thus, any path along 7 between 2 distinct points 0/7 (7 <9X is a generalized 
tangle jlffi/ , that is, a piecewise analytic path intersecting <9E transversally and points of 
non-analyticity are vertices of 7, equivalently, intersections of tangles. The subspace of 
H a f having the tangle property will be called Htangie- 

The tangle property is a rather natural assumption about states £ because curves of 
non-zero parameter length running in <9£ themselves between vertices automatically have 
infinite length with respect to 5 a b- This is not the case for curves between vertices which 
just approach <9S but otherwise lie in the interiour. 

Under this additional assumption, we have in the limit e — > 0, putting ( |3.13| ), ( |3.16| ) and 
the remaining pre-f actor of —2/k from ( |3.10| ) together 



E ADM (S)f, = -2h 2 K £ ix^/ 7 . (3.17) 

i>ed£UV(7) Vv 

But recognizing — iX v = J v as the total angular momentum of / 7 at v we finally find 

E ADM f J = 2h 2 K Y fW' ( 3 - 18 ) 
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Expression Q3.18| ) defines a self-consistent family of operators {Eadm^} of operators which 
can be extended to infinite graphs and thus to TCtangie provided that the number of punc- 
tures of 7 with <9£ is finite. This defines the domain of Eadm- 
A number of remarks are in order : 

• The operator J* is the infinitesimal generator of gauge transformations at v and 
therefore commutes with V v because the volume operator is gauge invariant proving 
that (|3.18|) is densely defined on the restricted state space. If we would not require 



the tangle property then while the volume integral in ( 3.10 ) gives essentially X v 
the surface integral does not, rather it gives X' v = X v — S e (o)=u,ecd£ which does 
not commute with V v and so may map / 7 into a linear combination of states some 
of which may acquire zero volume and so (ft.!8|) would blow up. Thus, the tangle 
property is sufficient for Eadm to be a densely defined operator. Although we do 
not have a proof, it is almost granted that X' v f 1 will contain zero volume eigenstates 
in its expansion and so the tangle property would also be necessary ! 



As a striking bonus of the tangle property we easily see that ( |3.18| ) indeed defines 
a positive-semidefinite self-adjoint operator on the Hilbert space TLtangie '■ First of 
all, since the volume operator and the X l v are defined in terms of the operators 
X\ it follows that the family of operators ( |3.18| ) is consistently defined because the 



X\ are. Next, since [A*,t^,] = we see that we may order ( [3.18D symmetrically 
involving only terms of the form J*i J* and so ( |3.18 ) defines a symmetric operator 



because it leaves the graph 7 invariant J|, [|. Finally, the Laplacian — A^ := J l v J l v 
has non-negative discrete spectrum j v (j v + 1) where j v is the spin of the contractor 
of the generalized (non-gauge- invariant) spin-network states into which / 7 can be 
decomposed. Moreover, since — A^ and V v commute, they can be simultanously 
diagonalized. Thus, if A„ is the eigenvalue of the volume operator then the simul- 
tanous eigenstate / 7 is an eigenstate of the ADM energy with eigenvalue J2v 

This provides an explicit diagonalization of Eadm on TCtangie and demonstrates that 
it is a self-adjoint operator on TCtangie- Now, since states in TCtangie are not diffeo- 
morphism averaged at spatial infinity the part of the graph which is responsible 
for the spectrum of the ADM energy is untouched by the diffeomorphism group 
averaging. Next, the map 17 that makes a diffeomorphism invariant state a solution 
of the Hamiltonian constraint is a generalized projector at each vertex of the graph 
separately |12|, [13|, [14], |18|]. Therefore, it is the identity map at those vertices of 



7 that lie in <9E because for a gauge transformation iV — > the action of the 
constraint at those vertices is trivially zero. Thus the part of the state responsible 
for the spectrum and the adjointness relations of Eadm is unchanged by the map 17 
and are thus preserved when we go to the physical Hilbert space TC p h ys = ^TCtangie- 
Moreover, the spectrum is entirely discrete and non-negative, thus we have proved 
A Quantum Positivity of Energy Theorem ! 

Astonishingly, the proof of this theorem turned out to be surprisingly simple : the 
proof of the classical positivity theorem is much more complicated and uses the 
boundary conditions and the Einstein equations at various stages. Why did we not 
need (a quantum analogue of) these assumptions ? The answer is that we actually 
did use them : We used them in the definition of an asymptotically flat state, in 
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particular, that the volume of the state be non-vanishing and that it be a solution 
of the quantum Einstein equations. Since the energy functional is really given by 
E(N) = H(N) + Eadai{N) where the Hamiltonian constraint only vanishes on 
a solution, likewise the ADM energy operator Eadm only represents energy if we 
apply it to a solution \1/ of H'(N)ty = 0. Here we have written the dual H'(N) on 
$' of H(N) because, as stated above, solutions \& lie actually in $' ||, [l4| . Finally, 



we used a regularization of the operator consisting 1) in the restriction of the space 
of states to functions of classical connections (this means here that they are smooth 
and decay at infinity as 1/r 2 ) and then 2) in the extension of the expression for the 
operator obtained to all of A. In the derivation of the operator we made crucial use 
of the fact that a classical connection decays at infinity. 

Now a subtle issue is the following : by definition a solution \1> satisfies H'(N)^f = 
for any N which vanishes at <9X. But how about lapse functions N that approach 
a constant value at dH ? It is now not a consequence of the formalism any longer 
that H'(N)^ = should hold, very much like in the case of the Gauss constraint. 
The only guideline of what to do is the classical theory and there it is indeed 
true that on classical solutions (A , E ) to the Einstein equations E(N) just equals 
Eadm(N)a=a ,e=e so that we will require that H'(N)ty = even for asymptotically 
constant lapse. 

Now, by definition the operator fj acts like the identity operator at <9£. Therefore we 
conclude that any physical state has the property that it is annihilated by H'(N) 
even for N = const, and moreover it is a linear combination of eigenstates with 
non-zero eigenvalue of the volume operator V v for each vertex of the graph 7 on 
which £ depends (with r/£ = such that v G <9£. One might suspect that the 
number of states that satisfy this condition is rather tiny but the contrary is the 
case : the volume operator has the particular property that it does not change the 
graph or the labelling of that graph with spin quantum numbers. Now there exist 
an infinite number of states which are annihilated by H'(N) just because the graph 
or its labellings is of a particular type (see £L4| where such states were labelled 
by "spin- webs", more precisely, "sources" of spin- webs) and so one can construct 
eigenstates of the volume operator of such states while they are still annihilated by 
H'(N). As an aside, this might shed some light on the issue of how to interprete 
those special solutions of the Hamiltonian constraint. 
These remarks are sufficient to show that then E' ADM (N) = Eadm{N). 

The fact that the ADM energy operator is essentially diagonal on spin-network 
states can be interpreted as saying that the spin-network representation is the non- 
linear Fock representation of quantum general relativity. Namely, if we compare 
the spin that the edges of a spin-network carry with the occupation number of 
momentum modes of, say, the free Maxwell field then we may interprete this spin 
essentially as the occupation number of a gravitational mode which is not labelled 
by a momentum but by an edge. 

Notice that each vertex of 7 which lies in <9£ must be at least three-valent (not four- 
valent because the function does not need to be gauge invariant at infinity) in order 
that it can have non-vanishing volume. Also notice that the gravitational energy is 
quantized in quanta of the Planck mass : the eigenvalues of the volume operator 
are multipla of £p so the eigenvalues of energy, according to (|3.18| ) are multipla of 
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k = m p . 



• For a classical connection every function at <9E is gauge invariant because it decays 
like 1/r 2 which results in a holonomy of order exp(i/r) — > 1, i.e., a trivial holonomy. 
In quantum theory this is lifted by the distributional nature of a connection, smooth 
connections are assigned zero volume by the Ashtekar Lewandowski measure in the 
space of distributional connections and are unimportant. 

• Gauge invariant states at <9E correspond to vanishing energy eigenvalue. Thus, en- 
ergy seems to sit at non-gauge invariant vertices. We may interprete this observation 
as follows : The gravitational energy in a state labelled by an (infinite) graph 7 is 
concentrated at the vertices of 7, and energy flows from vertex to vertex along the 
edges of 7 in quantized packages labelled by the spin of those edges. Non-zero en- 
ergy at a vertex corresponds to lack of gauge invariance at this vertex meaning that 
the spins that flow out or into a vertex do not add up to zero. Now in the interiour 
of E the Quantum Gauss Constraint requires that all spins add up to zero. We 
interprete this as the connection dynamics version of the geometrodynamics result 
that there is no energy location in general relativity in the interior of E, gravita- 
tional energy can be gauged away locally, it is pure gauge. However, while it can be 
pushed around at one's will, one cannot entirely delete it, one can push it all the 
way to spatial infinity where it eventually shows up in the form of a non-zero net 
spin flow at the vertices of 7 at <9E. 

• The fact that (representations of) the SU (2) gauge group of general relativity should 
play an essential role for the energy is very unexpected from a geometrodynamics 
point of view where one never even talks about the SU(2) gauge freedom. Even the 
classical ADM expression J as dS a (q a b,b — Qbb,a) is manifestly gauge invariant, so how 
did the SU{2) gauge group enter the stage?] ? The answer is the following : Notice 
that while the classical ADM expression is manifestly gauge invariant, it is not at 
all covariant : The derivatives that appear in J as dS a {q a b^ — Qbb,a) are not covariant 
derivatives (in fact if they were covariant derivatives then the energy would vanish 
identically). This does not need to concern us because at <9E the diffeomorphisms 
that underly the diffeomorphism constraint must die off. Now as we showed in ( ^.2| ), 
when writing J as dS a (q a bfi — <?w>,a) in terms of Ef this lack of three-diffeomorphism 
covariance, by means of the boundary conditions, gets translated into non-gauge- 
invariance while installing covariance because in ( |3.2| ) a correction term drops out 
due to the boundary conditions which explains why one and same function can 
be written in a manifestly gauge invariant or diffeomorphism covariant way but 



not both. Indeed, the expression —2 j QS dS a EfdbE\ / det(g) is manifestly three- 

diffeomorphism covariant but it fails to be gauge-invariant. This is not unexpected 
: after all the triad formulation reduces the local diffeomorphism gauge freedom to 
local rotation freedom. In conclusion, the fact that states with non-zero energy are 
not gauge invariant is in fact very natural. 

2 This question actually arises already in connection with the spectrum of the geometrical operators 
volume, area and length [ p7[ p8| , [2S| , [§0[ |l6| ], however, since irreducible representations also carry gauge 
invariant information the fact that these operators have a spectrum which is determined by spin-quantum 
numbers of gauge-invariant states is maybe not that surprising. What is surprising for the energy operator 
is that it is the spins of non-gauge invariant states which determines the spectrum. 
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Another function for which SU(2) gauge transformations and diffeomorphisms get 
mixed up is the classical vector constraint V a = tr (F a bE b ). Strictly speaking this 
constraint function does not generate diffeomorphisms but only on gauge invariant 
functions. 



Notice that although —2 / 9E dS a EfdbE b / y det(g) is not gauge invariant the quantum 
expression (|3.18|) in fact is. The reason why that is possible lies in the structure of 
quantum theory : in the classical theory we only have functions on phase space. In 
quantum theory those functions get translated into operators on a Hilbert space but 
values of those functions really correspond to expectation values. Thus non-gauge 
invariant functions correspond to expectation values of either a gauge-invariant op- 
erator in a non-gauge-invariant state or vice versa. The quantization (|3.18|) picks 
the latter possibility. 

In principle we have now solved the "problem of time" : Since we have a true 
Hamiltonian we can introduce the Schroedinger time parameter t and our state 
vectors \I/ £ <3>', being distributions which are invariant under asymptotically identity 
diffeomorphisms are supposed to satisfy the non- stationary Schroedinger equation 

Eadm(N = = -ihdt* . (3.19) 

Notice that, as we showed above, the ADM energy operator is its own dual so that 



( |3.19| ) makes sense. 



In the next section the operator (|3.18|) will be shown to commute with all quantum 



constraints of general relativity and it is therefore a Strong Quantum Dirac Observ- 
able for Quantum Gravity in the strict sense of the word. This is not unexpected 
because it is built purely from momentum operators. 



4 The Poincare Algebra 

We now wish to define the rest of the quantum generators of the little group of the 
asymptotic Poincare group and check whether their algebra is anomalous or not 0. This 
is enough to construct particle states since the irreducible unitary representations of the 
little group induce a unique unitary irreducible representation of the full Poincare group. 
So far we did not construct an operator corresponding to a boost generator which is more 
difficult to obtain than the ADM energy operator. 

First of all we must clarify on which space to represent the Poincare group, respec- 
tively its generators. To that end it is helpful to remember how the classical Poincare 



generators are realized as a subalgebra of the Poisson algebra [g(J, ^1 



3 By little group of the Poincare group we mean the group generated by the four-translations and the 
little subgroup of the connected component of the Lorentz Group. The latter, as is well-known, is the 
stabilator subgroup of the Lorentz group associated with a standard four vector s. In the massive case 
s 2 > the standard vector is associated with the spin of the particle in the rest frame and the covering 
group of the stabilator group is given by by SU (2). In the massless case the standard vector is associated 
with the helicity of the particle (spin in momentum direction) and the covering group of the stabilator 
group is given by by U(l), physically important representations being two-valued. Thus, the rotations at 
spatial infinity determine the unitary irreducible representation of the particle state in question. 
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Let H(N),V(N) be the Hamiltonian and diffeomorphism constraint functional respec- 
tively. Both functionals are integrals over S of local densities and both converge and 
are functionally differentiable only if the lapse and shift functions N, N vanish at <9E. 
In order to be able to describe the Poincare group corresponding to the asymptoti- 
cally constant or even diverging functions (x a is a cartesian frame at spatial infinity) 
N = a + x a x a , N a = a a + e abc 4>bX c where (a, a a ) is a four translation, <p a are rotation angles 
and x a are boost parameters, one proceeds as follows : let S be a bounded two-surface 
that is topologically a sphere and let B(S) be the (intersection of £ with the) closed ball 
such that dB(S) = S. For each S one defines E(N,S) := H(N,S) + E ADM (N,S) + 
B(N, S), P(N, S) := V(N, S) + Padm{N, S) where the parameter S means that volume 
integrals are restricted to B(S) only (a classical regularization of the divergent integrals) 
and the "counter-terms" Eadm(N, S), B(N, S), Padm(N, S) are the surface integrals de- 
fined in ^Of and correspond to ADM energy, boost and momentum. One can show that 
lim<^a S E(N, 5),lim s ^ es P(N, S) exist. Moreover, for each finite S, E(N, S),P(N,S) 
are functionally differentiable so that it is meaningful to compute the Possion brackets 

{E(M, S), E(M, S)} = P(q ab (MN, b - M )b N), S) 
{E(M,S),P(N,S)} = E(C^M,S) 
{P(M, S), P(N, S)} = P(JC^N, S) . (4.1) 

The crucial point is that one computes the Poisson brackets a) at finite S and b) on the 
full phase space and then takes the limit S — > <9£ or restricts to the constraint surface of 
the phase space (where H(N, S) = V(N, S) =0). Notice that the numerical value of, say, 
E(N, S) equals H(N, S) for a gauge transformation for which iV — > as S —>■ dH. On 
the other hand, on the constraint surface for a symmetry for which iV y4 as S — ► it 
equals a time translation or a boost respectively. A similar remark holds for P(N, S). One 
therefore interpretes (|4.1| ) as follows : if M, N are both pure gauge then the constraint 
algebra closes. If M is a symmetry and N pure gauge then energy (or boost generator) 
are gauge invariant. If M, N are both symmetry then time translations commute with 
each other, time translations and boosts give a spatial translation and a boost with a 
boost gives a rotation, in other words the symmetry algebra closes. 
In quantum theory we will therefore proceed as follows : 



Recall 0, 13, 14, 18] that the Hamiltonian constraint H(N) (for asymptotically vanishing 
N) is only well-defined on the subspace of $' corresponding to distributions on $ which 
are invariant under diffeomorphisms that approach identity at dH. Thus we can expect 
the symmetry algebra to hold only on such distributions as well. In fact, we will just 
choose ^ to be a solution to all constraints. 

Next, in view of the fact that even the classical symmtry algebra only holds provided 
one first computes Poisson brackets at finite S and then takes the limit, we will check 
the quantum algebra first by evaluating \1/ on E(Ns)fs for functions /j6$ cylindrical 
with respect to a graph which lies in the interior of B(S) (it may intersect S in such a 
way that the volume operator does not vanish at the intersection point for none of the 
eigenvectors into which fg maybe decomposed) and lapse functions N$ which grow at 
infinity like symmetries but which are supported in B(S) U S including S, and then to 
take the limit S — > <9£ (the support fills all of S as S — > <9£ in this process). 

We come to the definition of E(N), P(N). First we treat the spatial Euclidean group. 
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The unitary representation of the diffeomorphism group defined by U{ip)f 1 = which 
was for matters of solving the diffeomorphism constraint so far only defined for diffeo- 
morphisms that approach asymptotically the identity, can easily be extended to three- 
diffeomorphisms which correpond to asymptotic spatial translations or rotations. Instead 
of defining the generator P{N) though (which does not exist on 7i @) we content ourselves 
with the exponentiated version U((p(N)) where <f(N) is the diffeomorphism generated by 
the six parameter shift vector field N a = a a + e a b c <t> h x c for some cartesian frame x a possibly 
corrected by an asymptotically vanishing vector field corresponding to a gauge transfor- 
mation. It is trivial to check that 

^(^(iV))^^^))^^^))' 1 ^^^))- 1 = UMCjfN*)) (4.2) 

where £ denotes the Lie derivative so that there are no anomalies coming from the spatial 
Euclidean group. This expression was derived by applying it to any function f$ cylindrical 
with respect to a graph with support in B(S). 

We now turn to the time translations. As already mentioned we will not consider 
boosts in this paper so that Xa = in the four parameter family of lapse functions 
N = a + XaX a (modulo a correction which vanishes at <9£). Define the operator on 7i 

E(N) :=H(N) + E ADM (N) (4.3) 

where H(N) is the Lorentzian Hamiltonian constraint. Notice that E(N) just as the 



Hamiltonian constraint in fl2| , |13| , [14], [18] carries a certain prescription dependence which 
is removed by evaluating its dual on $d«//- We will not repeat these details here and 
refrain from indicating this prescription dependence in |4.3| , however, the prescription 
dependence has consequences for the commutator algebra that we will dicuss below in 
great detail. 

Let us verify the commutators between the time translations among themselves and 
between time taranslations and spatial translations and rotations. We have 

*([£(M), E(N)]f y ) = *([#(M), H(N)]f y ) + *([E ADM (M), E ADM {N)]f 7 ) 
+ V{{[E ADM {M),H{N)} + [H{M),E ADM {N)]}fJ . (4.4) 

The first term vanishes for the same reason as in |OJ, [H| FJ, |T^] although one needs one 
additional argument : the Hamiltonian constraint does not act at vertices that it creates. 
Therefore, it can be written as a double sum over vertices v, v' of 7 alone and each of 
these terms is of the form 

(M(v)N(v') - M(v')N(v))V([H v , Mv) H v „ - H^^M) 

where the notation means that H Vjl is a family of consistently defined operators each of 
which acting on cylindrical functions which depend on the graph 7 and 7(f) is a graph 
on which H va f 1 depends. This expression clearly is non-vanishing only if v ^ v' but 
then it can be shown that the operators H v> „ and H^,.. actually commute. Now still this 
does not show that the term above vanishes however, it can be shown that H v i i yt v \H va f 1 
and H V)1 ( v i\H v i a f 1 are related by a diffeomorphism ||13|| . Now in Jl3| that was enough 



to show that the commutator vanishes because we were dealing there only with vertices 
which do not intersect 5" as otherwise both lapse functions identically vanish for a pure 
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gauge transformation. Thus the diffeomorphism that relates the two terms above could 
be chosen to have support inside B(S) and \I/ is invariant under such diffeomorphisms. 
In the present context that does not need to be true. However, the crucial point is now 
that by the tangle property all edges of 7 that intersect S must intersect S transversally. 
Therefore the arcs that the Hamiltonian constraint attaches to 7 and whose position is the 
only thing by which the two above vectors differ lie inside B(S) and do not intersect S. 
Therefore, again the two vectors are related by a diffeomorphism which has support inside 
B(S), that is, they are related by a gauge transformation and therefore the commutator 
vanishes. 



We turn to the second term in (fO). Now we obtain a double sum over vertices of 7 



which lie in S and each term is of the form 

(M(v)N(v') - M(v')N(v))^([E v , ADM ,E vADM ]f,) 

which is significantly simpler than before because E v> adm does not alter the graph. Notice 
that the commutator makes sense because Eadm,v leaves the span of non-zero volume 
eigenvectors invariant. Now for v 7^ v' the commutator trivially vanishes, this time 
without employing diffeomorphism invariance of ^. 

Finally the last term in Q4.4T) is a double sum over vertices v, v' of 7, where v must lie 
in S, of the form 

(M(v)N(v') - M(v')N(v))^([H vt> „ E vADM ]f,) . (4.5) 



The fact that Eadm does not alter the graph was used to write (|4.5|) as a commutator 
without employing diffeomorphism invariance of Now it may happen that, although 
/ 7 is in the domain of E V) adm, that H va f 1 is not any longer in the domain and so (|4.5|) , 
for v = v', is in danger of being a meaningless product of something that blows up times 
zero while that cannot happen for v 7^ v'. However, since ^ is a solution we conclude first 
of all that ( |4.5| ) equals 

- (M(v)N(v') - M(v')N(v))[E vADM y](H v ,^f,) (4.6) 

and since $ is also in the domain of Eadm both E Vj adm^ an d H v > a f^ are well-defined 
elements of $' and $ respectively we conclude that in case v = v ' (|4.5|) indeed vanishes. 
On the other hand, the same argument as before shows that the commutator trivially 
vanishes for u^/ 

Let us now check the commutator between time translations and spatial translations 
and rotations tp. We have 

V(lU(<p)- l E(N)U(<p) - E(N)]f,) 

+ Mv(v)MU(<p- 1 )E A DMMv)U(<r)) ~ NivW&ADMM] . (4.7) 

«ev(7)ns 

Since Eadm does not change the graph on which a function depends we have identically 

U(tp~ l )EADM,<p(y)f(p(j) = EADM,vfj- 

Now, as explained in more detail in [0, the operator H(N) depends on a certain prescrip- 
tion of how to attach loops to graphs. Since in the interiour of B(S) there is no background 
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metric available, this prescription can only be topological in nature and therefore graphs 
differing by a diffeormorphism (p are assigned graphs by H(N) which are diffeomorphic 
by a diffeomorphism which may not coincide with (p. That is, in the interiour of B(S), 
H(N) is only covariant up to a diffeomorphism. On the other hand, since one has the 
fixed background metric 5 a b at S one can make H(N) precisely covariant at S, that is, the 
prescription satisfies tp\s = f\ s - Therefore, with this sense of covariance of H(N) it is true 

that C/(<^ -1 )H' v ,( v ) jV ,( 7 )/ v ,( 7 ) and H v>1 f 1 differ at most by a diffeomorphism with support in 
the interiour of B(S). 
In conclusion we obtain 

*([E(N), £%)]/ 7 ) = *(E(<p*N - A0/ 7 ) 

which is what we were looking for. 

We conclude that the little algebra of the Poincare algebra is faithfully implemented. 
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